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We report results of the first 3x3 "multi-angle" simulation of the evolution of neutrino flavor 
in the core collapse supernova environment. In particular, we follow neutrino flavor transformation 
in the neutronization neutrino burst of an O-Ne-Mg core collapse event. Though in a qualitative 
sense our results are consistent with those obtained in 3 x 3 single-angle simulations, at least in 
terms of neutrino mass hierarchy dependence, performing multi-angle calculations is found to reduce 
the adiabaticity of flavor evolution in the normal neutrino mass hierarchy, resulting in lower swap 
energies. Differences between single-angle and multi-angle results are largest for the normal neutrino 
mass hierarchy. Our simulations also show that current uncertainties in the measured mass-squared 
and mixing angle parameters translate into uncertainties in neutrino swap energies. Our results 
show that at low #13 it may be difficult to resolve the neutrino mass hierarchy using the O-Ne-Mg 
neutronization neutrino burst. 



PACS numbers: 14.60.Pq, 97.60.Bw 

I. INTRODUCTION 

Here we present results of the first multi-angle, self- 
consistent, 3-flavor simulation of neutrino flavor evolu- 
tion in the core collapse supernova environment. With 
an expected ultimate energy release of ~ 0.1 M as neu- 
trinos of all kinds, these supernovae have long been stud- 
ied as potential sources of neutrinos which could provide 
probes of physics beyond the reach or scope of conven- 
tional terrestrial experiments. Neutrino flavor oscilla- 
tions are a case in point. Though we already know a 
great deal about neutrino mass-squared differences and 
flavor mixing properties from experiments, there are still 
fundamental neutrino mixing physics unknowns, e.g., the 
neutrino mass hierarchy, mixing angle #13, and CP vio- 
lating phase S. Assessing how the known and unknown 
neutrino flavor oscillation physics affects the core collapse 
neutrino burst signature is a necessary step in under- 
standing supernovae, as well as in gleaning insight into 
fundamental neutrino physics. 

We have chosen to examine the neutrino emission 
from supernovae originating from stars in the mass range 
~ 8 — 10 Mq. Supernova progenitors in this mass range 
are expected to be relatively common, comprising 25% 
or more of core collapse supernova events (taking the 
number of stars per unit mass to be cx m~ 2 ' 35 , using 
the Salpeter initial mass function). These relatively light 
supernova progenitors develop O-Ne-Mg cores, which un- 
dergo gravitational collapse triggered by electron capture 
on Ne and Mg isotopes. This collapse leaves a thin enve- 
lope above the proto-neutron star remnant, in contrast to 



the envelopes of more massive stars that have undergone 
the gravitational collapse of an Fe core. 

Progenitors of O-Ne-Mg core-collapse supernovae are 
characterized by a steep matter density gradient above 
the core, which may have shed most of the hydrogen en- 
velope prior to core collapse due to pulsation and winds. 
These characteristics lead to two salient features of O- 
Ne-Mg core-collapse events. First, in contrast to the Fe 
core-collapse case, the bounce shock in models of these 
supernovae is not severely hindered by the pressure of 
material falling onto the core. Consequently, explosion by 
direct neutrino driven shock re-heateing is obtained [IH3 ■ 
Second, the steep matter gradient above the core allows 
neutrino self-coupling (neutral-current neutrino-neutrino 
forward exchange scattering) to engineer collective neu- 
trino flavor transformation [4-16J during the early phase 
of neutrino emission, including shock break-out and the 
attendant neutronization neutrino burst [T71 HH] . 

However, there is another consequence of the steep 
matter density gradient. A steep matter gradient ne- 
cessitates accounting for the physics of the interfer- 
ence between neutrino flavor mixing scales. Essentially, 
Mikheyev-Smirnov-Wolfenstein (MSW) resonances can 
occur very near to each other in physical and energy space 
when the matter density profile drops off rapidly with ra- 
dius. The near overlapping of resonances for the solar and 
atmospheric vacuum neutrino mass-squared differences, 
Am| = 7.6 x 10~ 5 cV 2 and Am 2 tm = 2.4 x 10~ 3 eV 2 , 
respectively, dictates that a full 3x3 flavor treatment 
must be employed. Because these resonances occur close 
to the core, and because in this case the number density 
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of neutrinos emitted from the core will fall much more 
slowly than the local matter density, the anisotropic na- 
ture of the supernova environment necessitates a multi- 
angle simulation that includes trajectory-dependent neu- 
trino self-coupling [T71 flSH^S] . 

Nonlinear neutrino flavor transformation in supernovae 
is a subject in its infancy and, by necessity, much of 
the discussion in the supernova community is focused 
on technical issues on how to calculate it. To that end 
this paper raises an important point, that a "standard" 
simplifying treatment of the problem, the widely used 
"single-angle" approximation, may not reproduce the re- 
sults of more sophisticated "multi-angle" calculations. 

Reference [T7] , which employs the single-angle approx- 
imation, found two flavor swaps in the normal mass hi- 
erarchy for the O-Ne-Mg core collapse neutronization 
burst, one for each mass-squared splitting, but only one 
flavor swap in the inverted mass hierarchy. The gen- 
eral features found in the results of Ref. [T7] seem to 
agree with the semi-analytical analysis in Ref. [H]. Ref- 
erence [27] has pointed out that in single-angle calcula- 
tions, under some circumstances, a single mass-squared 
splitting can give rise to two spectral swaps. However, 
Ref. [28] demonstrated recently that full 3x3 simulations, 
as opposed to a sequence of 2 x 2 level crossings, may be 
necessary to understand the formation or suppression of 
these multiple flavor swaps. 

As in Ref. [T7], to better understand multi-angle ef- 
fects in 3 x 3 neutrino flavor oscillation scenarios and to 
enable direct comparison between single-angle and multi- 
angle caculations, we have chosen to simulate the epoch 
of the neutronization neutrino burst. This corresponds 
to an epoch only some ~ 10 ms after core bounce, when 
the shock wave propagates through the neutrino sphere. 
The total neutrino luminosity at this epoch is comprised 
predominantly, but not exclusively, of electron type neu- 
trinos (y Q ) left over from core collapse-driven neutroniza- 
tion. 

In section II we discuss issues associated with 3- 
neutrino flavor oscillations in vacuum and in medium. In 
section III we discuss our particular numerical method- 
ology for simulating multi-angle, 3-neutrino flavor devel- 
opment in supernovae. In section IV we discuss results 
for the neutronization burst in O-Ne-Mg core collapse 
events. Conclusions and prospects for future work are 
discussed in section V. 



II. SUPERNOVA NEUTRINO FLAVOR 
TRANSFORMATION 

Neutrino oscillations in vacuum arise because the weak 
interaction (flavor) eigenstates for these particles \v a ) 
are not coincident with their energy (mass) eigenstates 
\vi). Here a = e,/i,r and i — 1,2,3 refer to the flavor 
states and mass eigenvalues m u respectively. The neu- 
trino mass-squared differences are Amjj = mf — m|. In 
vacuum the energy eigenstates are related to the flavor 



eigenstates by the Maki-Nakagawa-Sakata (MNS) ma- 
trix U: \i/ a ) = Y^i^aiWi)'' wri ere U a i are the elements 
of the unitary transformation matrix. This transforma- 
tion has four free parameters: three mixing angles and a 
CP-violating phase 5. Solar, atmospheric, reactor, and 
accelerator neutrino experiments have measured two of 
these parameters: # 12 ~ 0.59, 623 ~ 7r/4. The best cur- 
rent experimental constraints show that #13 < 0.2 (2a 
limit) [2H3- The CP- violating phase S remains uncon- 
strained. 

The ordering (hierarchy) of the neutrino mass-squared 
differences also remains undetermined. There are two 
possible configurations for a set of 3 mass states with 
two mass-squared splittings. The "normal" mass hierar- 
chy has the solar neutrino mass-squared split below the 
atmospheric split. In this case the vacuum mass states 
are ordered from lowest to highest as mi, ni2, and 7713, 
with Am| = ml - m\ and Am^ = m§ - m|. By con- 
trast, the "inverted" mass hierarchy is where the solar 
neutrino mass-squared split lies above the atmospheric 
split. In this scheme the vacuum mass states are or- 
dered from lowest to highest as m 3 , mi, and m 2 , with 
Am.Q = m| — m\ and Am^ tm = m\ — m|. 

Neutrino oscillations in medium can differ significantly 
from the vacuum case. For the supernova environment 
above the proto-neutron star the neutrinos are nearly all 
free streaming. In this regime we therefore can neglect 
inelastic processes and associated de-coherence. Conse- 
quently we follow only coherent, elastic neutrino interac- 
tions. In this limit we can make the mean-field, coherent 
forward-scattering approximation, where neutrino flavor 
evolution is governed by a Schrodinger-like equation of 
motion. For example, we can represent the flavor state 
of neutrino n by \ip Vin ). (The relationship between neu- 
trino and antineutrino flavor states and the correspond- 
ing flavor isospin spinors is discussed in Ref. [21].) The 
evolution of this flavor state in the mean field coherent 
limit is then 

t dt = HWv,n) (1) 

where t is an Affine parameter along neutrino n's world 
line and H is the appropriate flavor-changing Hamilto- 
nian along this trajectory: H = H vac + H ma ,t + H vu . 
Here H vac , H mat , and H vll are the vacuum, neutrino- 
electron/positron charged current forward exchange scat- 
tering, and neutrino-neutrino neutral current forward ex- 
change scattering (neutrino self-coupling) contributions, 
respectively, to the overall Hamiltonian. 

There are two different mass-squared difference scales 
as discussed above, the atmospheric Am^ tm and solar 
Am.Q splittings. In what follows we refer to neutrino 
mixing at a particular point as being on the "Am Q scale" 
or "Am^ tm scale". By this we mean that the neutrino 
flavor transformation at this point is taking place mostly 
through v c ^ Vfi T mixing in that part of the unitary 
transformation corresponding to Am| or Am^ tm , respec- 
tively. This terminology is a holdover from the standard 
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adiabatic MSW case, where Am@/2E U and Am 2 atnl /2E V 
essentially pick out density regions and neutrino ener- 
gies E v where mixing is large (i.e., MSW resonance at 
Am? /2E U ~ v^Gf^-o, where n c is the net electron num- 
ber density). 

The coherent neutrino-neutrino forward exchange scat- 
tering Hamiltonian H vv produces vexing nonlinear cou- 
pling of flavor histories for neutrinos on intersecting tra- 
jectories. This is the pivotal complication encountered 
when attempting to calculate the evolution of the su- 
pernova neutrino flavor field. Note that H vv gives rise to 
both flavor-diagonal and off-diagonal potentials. In turn, 
each of these potentials is neutrino intersection angle- 
dependent, reflecting the V — A structure of the underly- 
ing current-current weak interaction Hamiltonian in the 
low momentum transfer limit. For neutrinos in state I 
on a trajectory with unit tangent vector fc;, the neutrino 
self-coupling Hamiltonian is given by a sum over neutri- 
nos and antineutrinos m: 

H vu j = V2G F ^2 (l - k ■ km) n„ >m |Vw) (Vv, m | 

m 

- V2G F ^ f 1 - ' ^w) n e,"L \lpp,m) (i>9,m\ ( 2 ) 

m 

where k m is the unit trajectory tangent vector for neu- 
trino or antineutrino m. Here n„ lT7l is the local number 
density of neutrinos in state m. In the supernova envi- 
ronment, the effect of neutrino self-coupling may be to 
engender collective neutrino flavor transformation, the 
general features of which are reviewed in Ref. |26j . 

At a particular epoch, and with a specified matter den- 
sity profile and specified initial neutrino fluxes and en- 
ergy spectra, numerical calculation of the neutrino flavor 
field above the proto-neutron star would require a self- 
consistent solution of Eq. [I] for neutrinos on all trajecto- 
ries. This must include a prescription for treating H vv in 
Eq. [2j which couples flavor evolution on intersecting neu- 
trino trajectories. This can be done by adopting the so- 
called "single-angle" approximation [10] , where the flavor 
evolution along a specified neutrino trajectory {e.g., the 
radial trajectory) is taken to apply along all other tra- 
jectories at corresponding values of the Affine coordinate 
on those world lines. Alternatively, a "multi- angle" treat- 
ment can be employed where no approximations are used 
in the self-consistent evaluation of the Hamiltonian in Eq. 
[2] Though especially difficult to carry out computation- 
ally in a full 3-neutrino mixing scheme, some supernova 
scenarios may require such a treatment. To this end, we 
have developed numerical techniques to implement this 
approach. 

III. METHODOLOGY 

The numerical codes used for simulating neutrino fla- 
vor evolution in the calculations reported in this pa- 
per are the FLAT code and the BULB code. These 



codes, and related schemes to solve for the flavor evo- 
lution of core collapse supernova neutrinos, are discussed 
in Ref. [25]. 

In order to parallelize the nonlincarly-couplcd differen- 
tial equations which describe neutrino flavor evolution, 
BULB employs a specific geometric representation of the 
region above the neutrino sphere. In this representa- 
tion all neutrinos are assumed to be emitted from a hard 
spherical shell, and propagate through a one dimensional, 
spherically symmetric distribution of matter. Of course, 
spherical symmetry implies that no aspect of any coher- 
ent neutrino forward scattering potential depends on po- 
lar or azimuthal coordinate. This high degree of symme- 
try allows the neutrino emission to be broken down and 
grouped by species, energy, and emission angle. Here 
we define the emission angle, to be the angle be- 
tween the neutrino direction and the vector normal to 
the surface of the neutrino sphere at the neutrino emis- 
sion point. Our choice of parallelization is to assign all 
of the neutrino species and energies for a single emission 
angle bin to a single core (multiple angle bins can be 
assigned to a single core if desired, though this is less 
efficacious). This allows for efficient, fine grained paral- 
lelization. This feature is critical, as the emission angle 
dimension in the simulation may require high resolution 
in order to achieve convergence of the whole calculation. 

To initialize the simulation, neutrinos are allocated to 
each energy-angle bin according to the species-specific lu- 
minosity and neutrino energy spectral type (blackbody or 
"pinched"). From there, BULB employs Heun's method, 
which is a second order predictive-corrective algorithm, 
to compute the flavor evolution of the neutrino states. 
To minimize inter-node traffic, each process sums up the 
neutrino density matrix elements for the energy-angle 
bins it is responsible for, and sends this reduced data 
back to a central process. This central process gath- 
ers all of the neutrino density matrix data and com- 
putes a unique forward scattering Hamiltonian Hq, in- 
cluding iJvac, H mat , and H vv , for each energy-angle bin. 
These Hamiltonians are then redistributed to the ap- 
propriate processes so that each process can evaluate 
\ip)At = CX P {—iH At) |^i)o, where \ip} is the initial wave 
function, and \ip)At is the resultant wave function after a 
single step in Affine parameter At along a given neutrino 
trajectory. Once that is completed, the central process 
again collects all of the neutrino energy-angle states \4>)At 
to compute a new set of Hamiltonians H& t appropriate 
for the end point of the step. These are then sent back 
out to the individual processes, which now use the av- 
erage of Hq and H\ t to compute the evolved neutrino 
state: 

Mfiaai = \ [exp (-vff At) + exp (-iH At At)} |V>)o- (3) 

In order to check for convergence, a second round of 
computations are made using the Heun method to evolve 
the state \ij))o, using the above algorithm twice, with 
a step size At/2. Every process checks that both final 
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states for 1-0) final agree to within a predefined error tol- 
erance, usually chosen to be 1 part in 10 s . If the final 
states for single energy-angle bins do not converge, the 
step size is halved and the process is started over. This 
guarantees that the entire set of neutrino flavor states 
are tracked consistently. Periodically, At is increased by 
a factor of 2 to make sure that the simulation is operating 
with an efficiently chosen time step. 

Convergence of the overall calculation is checked by 
comparing results with different error tolerances and 
differing numbers of energy and angle bins. Gener- 
ally, large numbers of energy and angle bins, typically 
~ (1000 angle bins) x (400 energy bins), are required to 
achieve good results. 

Wc have validated our simulations by performing them 
with two numerically distinct codes, FLAT and BULB, 
and comparing our results |25j . Using the same set of 
initial conditions, both codes agree with each other at 
the level of 0.1% when comparing the final neutrino fla- 
vor states. FLAT works in analogous fashion to BULB, 
though these codes use different integration algorithms 
for handling the partial differential equations governing 
neutrino flavor evolution. This procedure has enabled 
the development of a full 3x3 multi-angle treatment for 
supernova neutrino flavor transformation. 



IV. RESULTS 
A. Overview 

The results of our multi-angle, 3x3 simulations in 
a particular case are shown in Fig. [T] for the inverted 
neutrino mass hierarchy. We have chosen the values 
013 = 0.1, 0i2 = 0.6, 023 = tt/4, Amli = 8.0 x 10~ 5 eV 2 , 
and Am 2 2 = —3.0 x 10 -3 eV 2 for this particular calcu- 
lation to facilitate comparison with the choice of mixing 
parameters in Refs. [T71 ITS] , 

All of our simulations of the neutronization neutrino 
pulse from an ONeMg core-collapse supernova begin at 
an initial radius of r = 900 km, where the matter density 
is still large enough that no flavor transformation has yet 
taken place. In order to compare our results directly with 
the results of Ref. [T7], we use the density profile from 
the pre-collapse calculations in Refs. [551 [35] in a post- 
bounce epoch. However, the epoch in which we perform 
the calculations is ~ 10 ms post-bounce, while the free- 
fall timescale at the radius (900 km) at which we begin 
our calculations is r ~ 100 ms, so there is not much time 
for non-homologous modifications to the density profile. 
Moreover, the material at this radius is likely falling more 
slowly than the free fall rate. In any case, the objective 
of our calculations is to compare single and multi-angle 
treatments of neutrino flavor transformation. Ultimately, 
neutrino flavor transformation simulations must be per- 
formed in a consistent supernova model. 

We model the neutrino emission as originating from a 
uniform sphere at a radius of i?„ = 60 km above the core 



of the proto-neutron star. R v is the calculated radius of 
the "neutrinosphere" , where the electron neutrino optical 
depth is equal to unity. It has been observed that the 
neutrinosphere radius is not a sharp surface, but instead 
is partially smeared out depending on neutrino energy as 
discussed in Ref. [31] . In the case of ONeMg core-collapse 
supernovae the matter density profile falls so swiftly with 
radius that this effect is small compared to the radius of 
the neutrinosphere itself and also the distance above the 
neutrino-sphere where significant flavor transformation 
takes place. 

In Subsection B, we present results of calculations 
with more recent and accurate values for neutrino mass- 
squared differences. Figure [T] presents the electron flavor 
neutrino survival probability Pv e v a , as a function of the 
cosine of the emission angle $r. It also shows the vac- 
uum mass basis angle-averaged neutrino energy spectra. 
All of these results are at a radius r = 5000 km. In the 
simulations used to generate this figure we have followed 
Ref. |17j , and we have approximated the flavor content of 
the neutronization burst for the O-Ne-Mg core-collapse 
as pure v e . 

Our simulations show that neutrino flavor evolution 
in the inverted mass hierarchy agrees with what was 
found in Ref. [T7] in broad brush and also seems to 
agree with the explanation offered in Ref. [TS]. A clear 
step- wise swap occurs between the vijv\ mass states at 
E v ~ HMeV. As with the single-angle results, there is 
no large-scale neutrino flavor transformation generated 
by mixing at the Am 2 tm scale. 

It has been noted that in the inverted neutrino mass 
hierarchy the neutrino flavor field is unstable, in anal- 
ogy to the way a pendulum balanced in the inverted po- 
sition would be [22l EH]- Flavor transformation in the 
other neutrino and anti-neutrino species is non-linearly 
coupled to mixing at the Am 2 tm scale through the neu- 
trino self-coupling potentials. This coupling provides the 
impetus that drives the Am 2 tm flavor isospins away from 
their unstable equilibrium. However, the lack of any neu- 
trino species besides v c in the calculation shown in Fig.[l] 
prevents any perturbation of the "inverted flavor pen- 
dulum," leaving it essentially balanced. The inclusion 
of other neutrino and anti-neutrino species, which have 
energy luminosities L ~ Q.lL Ve at this epoch, does not 
significantly affect these results and conclusions. Figure[2] 
shows that relatively small numbers (~ 10% admixture) 
of these neutrinos are not capable of destabilizing the 
Am 2 tm neutrino flavor mixing. 

Figure [3] shows the results of the multi-angle simula- 
tions in the normal mass hierarchy. We have used the 
same mixing parameters and initial conditions employed 
in the inverted neutrino mass hierarchy calculation, ex- 
cept that the sign of the atmospheric mass-squared dif- 
ference has been reversed, Am| 2 = +3.0 x 10 _3 eV 2 . 
Again, these calculations agree in broad brush with the 
results obtained for the single-angle simulations reported 
in Ref. |17j . In our multi-angle simulations the v^jv<x 
swap occurs at an energy of E v ss lOMeV. This swap 
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FIG. 1: Left panel: electron neutrino survival probability P Ve v s (color/shading key at top left) for the inverted mass hierarchy 
is shown as a function of cosine of emission angle, cosier, and neutrino energy, E in MeV, plotted at a radius of r = 5000km. 
Right: mass basis (key top right, inset) emission angle-averaged neutrino energy distribution functions versus neutrino energy, 
E. The dashed curve gives the initial v a emission angle-averaged energy spectrum. Movies of this simulation can be found at 
the URL in Ref. [30] . Each frame of the movie shows a representation of the neutrino survival probability in various different 
bases at a fixed radius above the core. Each successive frame is 1 km further out from the initial radius of r; n i t = 900 km out 
to the final radius r — 5000 km. 
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FIG. 2: Emission angle-averaged neutrino energy distribution 
functions versus neutrino energy plotted in the neutrino mass 
basis for the 3x3 multi-angle calculation of neutrino flavor 
evolution. Results shown at a radius of r = 5000km. In this 
simulation, a small (10%) admixture of all other species of 
neutrinos and anti-neutrinos are included. 



arises from the Am^ tm splitting. Our simulations show 
that the second swap, corresponding to v<xjv\, is at an en- 
ergy of E v s» 13.5 MeV. This swap arises from collective 
oscillations generated by mixing at the Ami scale. 

In the case of the multi-angle calculation, both of 
the flavor swaps observed in the normal neutrino mass 
hierarchy occur at lower energies than their counter- 



parts in the single-angle calculations. Comparing the 
different simulations we find that for the ^3/^2 swap 
£wj£ i_anEle ~ 10.0 MeV and E^gf-angic _ 12.5 MeV, 

swap swap 

and for the v 2 /v 1 swap £^ i-an8lB « 13.5 MeV and 

£ S w n a S p C_anglc ~ 15.0 MeV. 

Interestingly, the inverted neutrino mass hierarchy 
does not exhibit any change in the observed v<z/v\ swap 
energy between the multi-angle and the single-angle cal- 
culations. The spectra we observe for the multi-angle cal- 
culations in the inverted neutrino mass hierarchy agree 
with the framework developed by Ref. [18 . In that work, 
a single-angle toy model calculation showed that the re- 
sultant spectra of the neutronization neutrino burst for 
an O-Ne-Mg core-collapse supernova was formed by an 
initial phase of non-adiabatic, synchronous MSW flavor 
transformation, followed by collective neutrino oscilla- 
tions that form the swaps observed in the simulations. 

The normal neutrino mass hierarchy case is quite dif- 
ferent. The final spectra in Fig. [3] show that there are 
~ 10% fewer neutrinos that remain in mass state 3. This 
implies that the flavor evolution in the multi-angle cal- 
culation was less adiabatic than that in the single-angle 
calculation. 

To illustrate the importance of using full 3x3 flavor 
mixing in the case of an O-Ne-Mg core-collapse super- 
nova, we have performed a 2 x 2 multi-angle calculation 
for each mixing scale in the normal neutrino mass hierar- 
chy. Briefly, the swap energy of the Am^ tm mixing scale 
is -E SW ap ~ 10.0 MeV in our 2x2 multi-angle calculation, 
which corresponds closely with the swap energy found 
in the 3x3 multi-angle flavor mixing case. However, 
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FIG. 3: Left panel: electron neutrino survival probability Pv e v e (color/shading key at top left) for the normal mass hierarchy 
is shown as a function of cosine of emission angle, cosier, and neutrino energy, E in MeV, plotted at a radius of r = 5000km. 
Right: mass basis (key top right, inset) emission angle-averaged neutrino energy distribution functions versus neutrino energy, 
E. The dashed curve gives the initial v a emission angle-averaged energy spectrum. Movies of this simulation can be found at 
the URL in Ref. [30] . Each frame of the movie shows a representation of the neutrino survival probability in various different 
bases at a fixed radius above the core. Each successive frame is 1 km further out from the initial radius of r; n i t = 900 km out 
to the final radius r — 5000 km. 



the swap energy of the Aijiq mixing scale is found to 
be inswap ~ 11.0 MeV for the 2x2 multi-angle calcula- 
tion, which is lower than the corresponding vijv\ swap 
£ swap 13.5 MeV of the 3x3 multi-angle calculation. 

B. Sensitivity to Neutrino Mass-Squared 
Differences 

Since the publication of the original set of papers on 
the neutronization neutrino burst of an O-Ne-Mg core- 
collapse supernova [T71 [TS], the experimental limits on 
Am| and Am 2 tm have been refined, resting currently 
at the values of Anig = 7.6 x 1CT 5 eV 2 and Am 2 tm = 
2.4 x lCT 3 eV 2 [35j. We have conducted a set of 3 x 3 
multi-angle calculations to ascertain the effect the new 
experimental constraints would have on the flavor evo- 
lution history of the neutrinos released in the neutron- 
ization burst. Figure [4] shows the results of our calcula- 
tions for neutrino flavor transformation in the inverted 
neutrino mass hierarchy. There is no significant change 
observed in the final state of neutrinos at a radius of 
r = 5000 km relative to our earlier calculations with 
Am§! = 8.0 x 10~ 5 eV 2 . 

Figure [5] shows the results of our calculations for neu- 
trino flavor transformation in the normal neutrino mass 
hierarchy with the latest mass-squared values. With the 
new neutrino mixing parameters there is an increase in 
the observed swap energies for mixing at both the Attiq 
and Am 2 tm scales over the cases with the original mixing 
parameters adopted in |17j . The 1/3 jvi swap occurs at an 



energy of E v f=s 16.5 MeV, while the i>i /vi swap is pushed 
even higher to an energy of E v ~ 19.0 MeV, as opposed 
to E v « 10.0 MeV and E v « 13.5 MeV respectively. 

Because detection and characterization of supernova 
neutrinos is the ultimate aim of research in this field, 
we have included Fig. [6] and Fig. [7] to show predictions 
for the flux of electron neutrinos emitted by the neu- 
tronization neutrino burst of an O-Ne-Mg core collapse 
supernova. Figure [6] shows the predictions of both the 
single-angle and multi-angle calculations for the normal 
neutrino mass hierarchy, while Fig. [7] shows the same for 
the inverted neutrino mass hierarchy. 



C. Variation of #13 

Of particular interest is the effect produced by a signif- 
icantly decreased #13 mixing angle. As has been found in 
previous work, reducing the value of #13 can affect sensi- 
tively the energy of the flavor swap at the Am 2 tm scale for 
the normal neutrino mass hierarchy |24j . We performed 
a new calculation that matched the initial conditions of 
our primary, v c only, simulation of the neutronization 
burst, but this time with a value of 6*13 = 1.0 x 10~ 3 , 
two orders of magnitude smaller than the #13 value in 
the original calculation. This new value of # 13 drives 
the neutrino background-enhanced flavor evolution at the 
Am 2 tm scale entirely non-adiabatic, and Fig. [8] shows 
that this produces a resultant neutrino energy spectrum 
that is essentially indistinguishable from the spectrum 
produced in our calculations for the inverted neutrino 
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FIG. 4: This calculation of the flavor evolution of neutrinos in the inverted mass hierarchy was conducted with values of 
' eV 2 and Am 2 tm = —2.4 x 10 _3 eV 2 . Left panel: electron neutrino survival probability P VcVc (color/shading 



Am|, = 7.6 x 10" 



key at top left) for the inverted mass hierarchy is shown as a function of cosine of emission angle, cosier, and neutrino energy, 
E in MeV, plotted at a radius of r = 5000 km. Right: mass basis (key top right, inset) emission angle-averaged neutrino 
energy distribution functions versus neutrino energy, E. The dashed curve gives the initial v e emission angle- averaged energy 
spectrum. 
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FIG. 5: This calculation of the flavor evolution of neutrinos in the normal mass hierarchy was conducted with values of 
Am§) = 7.6 x 10 -5 eV 2 and Am 2 tm = 2.4 x 10~ 3 eV 2 . Left panel: electron neutrino survival probability P VtVe (color/shading 
key at top left) for the normal mass hierarchy is shown as a function of cosine of emission angle, cos$r, and neutrino energy, 
E in MeV, plotted at a radius of r = 5000 km. Right: mass basis (key top right, inset) emission angle-averaged neutrino 
energy distribution functions versus neutrino energy, E. The dashed curve gives the initial v c emission angle- averaged energy 
spectrum. 



mass hierarchy. This may present complications in the 
future determination of the neutrino mass hierarchy and 
6*i3 using the neutronization pulse of an O-Ne-Mg core- 
collapse supernova. Planned earth-based experiments 
should be able to set an upper limit on the value of 
sin 2 2^i3 ~ 1.0 x 10 -2 [35]. The neutronization neu- 



trino pulse signal from an O-Ne-Mg core-collapse super- 
nova may have limited ability to discern a value of #13 
much below this threshold. On the other hand, the late- 
time supernova neutrino signal, which should be more 
or less generic for all core-collapse supernovae, has no 
such limitations. Detection of the sense of the swap in 
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FIG. 6: Emission angle-averaged electron neutrino flux <3?„ 
(key top right, inset) for the normal neutrino mass hierarchy 
is shown as a function of neutrino energy E in MeV. The 
dashed curve gives the initial v E emission angle-averaged neu- 
trino flux. The shaded region gives the predicted flux in a 
single-angle calculation, and the thick line shows the flux pre- 
dicted by the multi-angle calculation. These calculations of 
electron neutrino flux are done using Am 2 } = 7.6 x ltr 5 eV 2 
10" 3 eV 2 and 0i 3 = 0.1. 
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FIG. 7: Emission angle-averaged electron neutrino flux 
(key top right, inset) for the inverted neutrino mass hierar- 
chy is shown as a function of neutrino energy E in MeV. The 
dashed curve gives the initial v E emission angle-averaged neu- 
trino flux. The shaded region gives the predicted flux in a 
single-angle calculation, and the thick line shows the flux pre- 
dicted by the multi-angle calculation. These calculations of 
electron neutrino flux are done using Adiq = 7.6 x 10 -5 eV 2 
and Am 2 tm = -2.4 x 10~ 3 eV 2 and 8 1S = 0.1. 



FIG. 8: Emission angle- averaged neutrino energy distribution 
functions versus neutrino energy plotted in the neutrino mass 
basis for the 3x3 multi-angle calculation of neutrino flavor 
evolution. This calculation employs a significantly reduced 
value of #13 = 1.0 x 10 -3 . This mixing angle is associated with 
flavor transformation at the Am 2 tm scale. Results shown at 
a radius of r = 5000km. 



this case should provide an unambiguous determination 
of the neutrino mass hierarchy. Measurements like these 
could be complementary to laboratory based neutrino 
mass probes. For example, experiments such as the MA- 
JORANA project may be able to determine the neutrino 
mass hierarchy directly |37j . 



V. CONCLUSION 

We have reported the first results of a fully 3x3, 
multi-angle simulation of the flavor evolution of neutrinos 
emitted from a supernova event. This case study exam- 
ined the neutrino flavor oscillations in the neutronization 
pulse emitted from an O-Ne-Mg core-collapse supernova. 
In these events, the matter density above the core has 
such a steep gradient that both neutrino self-coupling 
and the overlap of neutrino MSW resonance regions must 
be accounted for to obtain accurate energy spectra of the 
neutrino signals. 

We have developed a suite of programs that are capable 
of computing the solution to the non-linear, Schrddinger- 
like equation of motion that governs the neutrino fla- 
vor evolution in the coherent regime. As a result of 
the nature of the O-Ne-Mg core-collapse supernova envi- 
ronment, our codes have been pushed to a new level of 
complexity, handling the evolution of all three neutrino 
flavors with the inclusion of a full treatment of the an- 
gular dependence of both the neutrino-neutrino neutral 
current forward exchange scattering and neutrino flavor 
evolution history. By approaching the solution to this 
problem with two independent sets of solution algorithms 
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(the BULB and FLAT codes), we have verified that the 
results of our calculations are entirely self-consistent. 

The features of the neutrino flavor transformation his- 
tory for the neutronization neutrino burst observed in our 
calculations are partially consistent with what has been 
found in single-angle calculations. In the normal neutrino 
mass hierarchy, for # 13 = 0.1, a pair of flavor swaps are 
observed, corresponding to the Attiq and Am 2 tm mixing 
scales. In contrast, for this $13, the inverted neutrino 
mass hierarchy shows only a single flavor swap, origi- 
nating from collective neutrino oscillations at the AttIq 
scale. Our studies show that when #13 is reduced to the 
value of 0i3 = 0.001, the double swap observed in the 
normal neutrino mass hierarchy becomes a single flavor 
swap generated by mixing at the Amji scale, identical to 
the observed neutrino flavor energy spectrum calculated 
for the inverted neutrino mass hierarchy. 

In the normal neutrino mass hierarchy in particular, it 
is critical that 3x3 flavor evolution be employed. More- 
over, we have also shown that our results for the neu- 
trino flavor swap energies retain their relevance through- 
out the evolution of the neutronization neutrino burst it- 
self, as the added luminosity of neutrino species besides 
v Cl which rise late in the duration of the burst, do not 
significantly affect the energy of the swaps. 

We have shown that the energy of the flavor swaps can 
be sensitive to the values of Am 2 at the ~ 10% level. This 
gives added importance to continuing laboratory efforts 
to measure Am 2 . These constraints will only improve 
with time. 

It has been a hope that the collective neutrino oscil- 
lation signatures, if detected, could reveal the neutrino 
mass hierarchy and/or #13. However, discerning the neu- 
trino mass hierarchy from an O-Ne-Mg core-collapse neu- 
tronization burst signal would depend on distinguishing 
two spectral swaps from one spectral swap. Moreover, 
this would have to be done at relatively low neutrino 
energies, and likely also require both charged and neu- 
tral current detection capabilities. It is not clear that all 
proposed future supernova neutrino detectors would be 
capable of this. Our simulations also show that the nor- 
mal mass hierarchy will give rise to only a single swap 
if #13 is sufficiently small. This could make resolving 
the hierarchy problematic in a detected O-Ne-Mg core- 
collapse neutronization burst, depending on the value of 



#13. However, this limitation is likely not present in the 
late time neutrino signals of any core-collapse event. 

From the unexpected behavior of the swap energies 
in the normal neutrino mass hierarchy, it is clear that 
(as has been seen many times before in computational 
physics) incorporation of more realistic description of 
the physical environment in our simulations has resulted 
in new phenomenology. There is an intimate interplay 
at work between the competing effects of neutrino self- 
coupling and matter driven flavor evolution taking place 
in the heart of every exploding star. While the single- 
angle approximation remains a well understood and use- 
ful tool in evaluating the flavor evolution of neutrinos 
emitted from supernovae with relatively dense envelopes, 
those types of events may not be as common as previously 
believed. Given recent advances in the understanding 
of the hydrodynamics that governs core-collapse super- 
novae, it is clear that any hope of understanding the neu- 
trino signal from stars with thin envelopes may require 
even more detailed treatment [38l - l43j . 

In a sense our work is a step backward. At least 
for the particular supernova model we consider here, we 
must conclude that single-angle treatments are inade- 
quate when neutrinos have the normal mass hierarchy. 
It is still mysterious why the single-angle approach is ef- 
ficacious in some regimes and inadequate in others. In 
any case, the work presented here may be a first step 
towards understanding this issue. 
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